Direct mode summation for the Casimir energy of a solid ball 
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The Casimir energy of a solid ball placed in an infinite medium is calculated by a direct frequency 
summation using the contour integration. It is assumed that the permittivity and permeability of 
the ball and medium satisfy the condition eifii = £2/^2- Upon deriving the general expression for 
the Casimir energy, a dilute compact ball is considered (ei — £2) 2 /(ei + £2) 2 <C 1. In this case the 
QO ■ calculations are carried out which are of the first order in £ 2 and take account of the five terms in 

■ the Debye expansion of the Bessel functions involved. The implication of the obtained results to the 

ON \ attempts of explaining the sonoluminescence via the Casimir effect is shortly discussed. 
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!>• ' I. INTRODUCTION 

(N ; 

£SJ , The Casimir energy, determined by the first quantum correction to the ground state of a quantum field system 
^ ■ with allowance for nontrivial boundary conditions, proves to be essential in many problems of the elementary particle 
theory, in quantum cosmology, and in physics of condensed matter. However, up to now there is no universal 
method for calculating the Casimir effect for arbitrary boundary conditions. This has been done only for simple field 
configurations of high symmetry: gap between two plates, sphere, cylinder, wedge and so on. The curvature of the 
boundary and account of the dielectric and magnetic properties of the medium lead to considerable complications. 
While the attractive force between two uncharged metal plates has been calculated by Casimir as far back as 1948 [Q, 
this effect for perfectly conducting spherical shell in vacuum was computed by Boyer only in 1968 [|J (see also the 
latter calculations If an infinitely thin spherical shell separates media with arbitrary dielectric (£1,62) and 

magnetic (/ixj/-^) characteristics, this problem is not solved till now |^-|Io[|. The main drawback here is the lack of a 
consistent method for removing the divergences. Besides an attempt to revive the quasiclassical model of an extended 
electron proposed by Casimir fll|| , interest in this problem was also initiated by investigations of the bag models in 
C"| , hadron physics |l2|-|l4| and recently by search for the mechanism of sonoluminescence . 

In this paper we calculate the Casimir energy of a solid ball by making use of the direct summation of eigenfrequencies 
of vacuum electromagnetic field by contour integration p6| , pT| . A definite advantage of this method, compared with 
the Green's function technique employed in Q |9|,|l8| , is its simplicity and visualization. We consider a compact ball 
placed in an infinite medium when eifii = £2^2- This condition enables one to treat the divergencies analogously to 
the case of a perfectly conducting spherical shell jq] . Upon deriving the general expression for the Casimir energy, we 
address ourselves to the case of a dilute ball £ 2 <C 1, £ = (e± — 62) /(si + £2). The calculations here are of the first 
order in £ 2 and take account of the five terms in the Debye expansion of the Bessel functions involved. In this way 
we attain some generalization and refinement of the results obtained in this problem earlier |l8f| . 

The layout of the paper is as follows. In Sect. II we derive a general expression for the Casimir energy of a solid 
ball in an infinite surrounding under the condition £i/Ui = £2^2 — c~ 2 , where c is an arbitrary constant not necessary 
equal to one (it is the light velocity in the medium), the mode-by-mode summation of eigenfrequencies being used. 
In Sect. Ill the Casimir energy of a dilute ball £ 2 1 is calculated. The implication of the obtained result to the 
Schwinger attempt to explain the sonoluminescence via the Casimir effect is also considered. In Conclusion (Sect. 
IV) the results of the paper are briefly discussed. Dispersive effects are ignored in our paper. 
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II. CASIMIR ENERGY OF A SOLID BALL UNDER THE CONDITION ei/*i = £2^2 



Let us consider the Casimir of a solid ball of radius a, consisting of a material which is characterized by permittivity 
£\ and permeability fi\. We assume that the ball is placed in an infinite medium with permittivity £2 and permeability 
H2- We also suppose that the conductivity of the ball material and its surroundings is equal to zero. 

In our consideration the main part will be played by equations determining the eigenfrequencies u> of the electromag- 
netic oscillations for this configuration |ll| . It is convenient to rewrite these equations in terms of the Riccati-Bessel 
functions 



§i(x) = xji(x), e t (x) = xh^l^x), 



(2.1) 



where ji{x) — ^ir/2x Ji+\/2{x) is the spherical Bessel function and h^\x) — n/2xH^ 1 ^ 2 (x) is the spherical Hankel 
function of the first kind. For the TE-modes the frequency equation reads 



A ; (auj) = y/eifi2 s' l (k 1 a)ei(k 2 a) - ^/e 2 ^i sz(/cia)ej(fc 2 a) = 0, 



(2.2) 



where ki — ^j£i[ii lu, i = 1,2 are the wave numbers inside and outside the ball, respectively; prime stands for the 
differentiation with respect to the argument {k\a or k2a) of the corresponding Riccati-Bessel function. The frequencies 
of the TM-modes are determined by 



A™ (aw) = ^/£2Mi sj(fcia)ez(fc 2 a) - ^£1^2 sz(fcia)ej(fc 2 a) = 0. 



Under mutual change £i «-> fit , 



(2.3) 
= 1,2 



The orbital quantum num ber I in (2.2) and (2.3) assumes the values 1,2, 
frequency equations (|]^) and ( [2.3 ) transform into each other. 

It is worth noting that the frequencies of the electromagnetic oscillations determined by Eq. ( |2.2j ) and (2.3) are 
the same inside and outside the ball. The physical reason for this is that photons do not perform work when passing 
through the boundary at r = a. This is in contrast to the case of perfectly conducting spherical shell in vacuum ||, 
where eigenfrequencies inside the shell and outside it are determined by different equations |ll| . 

As usual we define the Casimir energy by the formula 



E 



5D 



■ LO, 



(2.4) 



where uj s are the roots of Eqs. (2.2) and (2.3) and u> s are the same roots under condition a — ► 00. Here s is a collective 
index that stands for a complete set of indices for the roots of Eqs. (2.2) and (|2.3|). Denoting the roots of Eqs. (2.2) 



and (2.3) by ct>£j and u))jV , respectively, we can cast Eq. (2.4) in the explicit form 



,(2) 



2 00 I 00 

^= EEEE( 

In=l 



,(«) 



a — 1 l — l m- 



1=1 



(2.5) 



where the notation 



z 00 



,(<*) 

nl 



-.(«) 

■'nl 



a— 1 n—1 



(2.6) 



is introduced. Here we have taken into account that the eigenfrequencies u>ff do not depend on the azimuthal 
quantum number to. For partial energies Ei we use representation in terms of the contour integral provided by the 
Cauchy theorem pO] 



Ei = 



1 + 1/2 
2-ixi 



c 



cl Aj E (az)Aj M (az) 
dz A™(oo)A™ (oo )- 



(2.7) 



where the contour C surrounds, counterclockwise, the roots of the frequency equations in the right half-plane. Location 
of the roots of Eqs. (2.2) and (2.3) enables one to deform the contour C into a segment of the imaginary axis (— iA, zA) 
and a semicircle of radius A in right half-plane. At a given value of A a finite number of the roots of frequency equations 
is taken into account. Thus A plays the role of a regularization parameter for the initial sum in Eq. ( |2.6| ) which sh ould 
be subsequently removed to infinity. In this limit the contribution of the semicircle of radius A into integral (2.7) 



2 



vanishes. From physical considerations it is clear that multiplier z in (2.7) is understood to be the lim^^o \/ z 2 + H 2 , 
where \x is the photon mass. Therefore in the integral along the segment (— iA,£A) we can integrate once by parts, 
the nonintegral terms being canceled. As a result Eq. (EJ) acquires the form 



1 + 1/2 [ Aj%y)A™( ly ) 
1 ira J V ATE(ioo)A™(joo)' 



(2.8) 



Now we need the modihed Riccati-Bessel functions 



/ 7TX I *2iX 

si(x) = J —I v (x), ei(x) = J— Kv(x), v = 1 + 1/2, 



(2.9) 



where I v (x) and K v {x) are the modified Bessel functions plfl . With allowance for the asymptotics of si(x) and ei(x) 
at x — > oo and fixed I 



si{x) ~ ~e x , 
ej(a;) ~ e _x 



(2-10) 
(2.H) 



equation (E^q) can be rewritten as 



i + 1/2 



dy In 



4 e -2(?i-(/2) 



(v^2 + V^T) 2 

u 

x [Vei£2MiA*2 ((sK<7iM<72)) 2 + M<7i)e!(g 2 )) 2 ) 
-(eiM2 + £2^i)si(gi)si(gi)ei(g2)eKq , 2)] 



(2-12) 



where qi — Je^JTi y, i = 1,2. We shall use this general equation in the next Section but here we address ourselves 
to the special case when the condition 



eiMi = £2^2 = c 



(2-13) 



is fulfilled. Here c is an arbitrary positive constant (the light velocity in medium). Physical implications of this 
condition at c = 1 can be found in J22| . Now Eq. ( 2.12| ) is simplified considerably 



El = <l±H21 ldyb > 



£ + £ 



(sKyMytf + is^eKyjf 



-(£ + £ 1 )si(y)s' l (y)ei{y)e' l {y) 



(2.14) 



where e — £\/£2- The argument of the logarithm in (2.14) can be transformed, if the following two equalities for the 
functions si(y) and ei(y) 



s'i(v) e l(y) - si(y)e[(y) = 1, 
s'i(y)ei(y) + siiyy^y) = (s ; (j/)e,(j/))'. 



(2.15) 
(2.16) 



are taken into account. It gives 



c(l+ 1/2) 



^ln{l-£ 2 [(siiyMy))'] 2 } , 



(2.17) 



where 



t M2 - Ml £ 1 _ £ 2 -2 ■ i o 

4 = ; = ; j e l \i l = c ,2 = 1,2. 

[12 + Mi £i + £2 



(2.18) 
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Thus , for a ball with a vacuum on the outside, £ = (1 — /x)/(l + //) = (e — l)/(e + 1) and c = 1. The expression 
( [2. 17 ) agrees with the results obtained in J§|,^3|], if one performs a partial integration of the expression for E given in 
thes e refe rences and puts the cutoff parameter 8 equal to zero. If e = or oo and c = 1 then, as one could expect, 
Eq. ( 2.17 ) turns into the analogous expression for the perfectly conducting spherical shell in vacuum H^]. 

We remark that in a prev ious paper (Tsl (1982)] an expression for the Casimir energy was calculated that seemingly 
is in conflict with Eq. ( 2.17 ). Namely, Eq. (2.42) in that paper corresponds to the following expression for Ei, assuming 
fii = fi, fi2 = 1 as above, and putting the cutoff parameter equal to zero, 



Here A; = (siei)', and Di, D\ are denned by 



DiDi dx 



Di(x) = i±si(x)e[{x) - s[(x)ei(x), 
Di(x) = fis'^eiix) - si(x)e[(x). 



(2.19) 



(2.20) 
(2.21) 



It turns out, however, that these two equations (2.17) and ( |2.19| ) are in agreement. To show the equivalence is not quite 
trivial, but follows after some algebra taki ng int o account the derivatives of the logarithms and the Wronskian (2.15). 
we omit the details here. The expression (2.17) is the most convenient form to work with, and we will use it in the 
following where we t urn to the special case of a dilute medium satisfying the condition e/i = 1 . We emphasize that 
the expression (2.17) is general, making no restriction at all on the magnitude of 



III. CASIMIR ENERGY OF A DILUTE COMPACT BALL 



Now we address ourselves to the consideration of the Casimir energy of a dilute compact ball when 

\t\ « I- (3.1) 

That means, to the lowest order in £, 



ln(l - £ 2 A 2 ) * ~£ 2 A 



/ > 



(3.2) 



which reflects a general property of all Casimir calculations in dilute media: the lowest order correction for all physical 
quantities is proportional to the s quare of the susceptibility (electric or magnetic). We shall henceforth work only to 
the second order in £. From (2.17) we then get for the Casimir energy 



E : 



i 

7ra 



2 00 i-oo 

- vl dx Xf(x). 



(3.3) 



For simplicity we have put here the constant c in condition (2.13) to be equal 1. 

We now invoke the following useful expansion, which was worked out by one of us some time ago [24| 



2v 



1- ^(2-27i 2 + 60i 4 -35i 6 ) 



128^ 



4(108 r - 3615 r + 21420 r - 47250 r 



+ 44352t 10 - 15015£ 12 ) + 0(l/i 



(3.4) 



Here, t(z) = (1 + z 2 ) 



2\-l/2 



functions [El|. From (3.4) we calculate 



xjv. This expression is based upon the Debye expansions for the modified Bessel 



A? 



Au 2 



1 



Av 2 



(2-27£ 2 + 60i 4 -35i 6 ) 



1 



■(1 - 54 i 2 + 1146 1 4 - 6200 i 6 + 13185 1 8 



16i> 4 

12138 1 10 + 4060 1 12 ) +0{l/v & ) 



(3.5) 
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which can now be inserted into Eq. (3.3). From the integral representation of the beta function [|TJ ~B(q,p) 
T(q)T(p)/T(q + p) we derive the formula 



p-l 



(3.6) 



which is quite useful in the present context. After some calculation we obtain 



E : 



3^ 
64a 



E" 

.1=1 



2-^ ,A + 1 R3R4 2^ ,A +U [ ,,< 



128 ^ ^ 2 ' 16384 ^ i/ 4 
z=i i=i 



(3.7) 



Here, the first sum will be regularized by means of the Riemann zeta function pq j, which turns out to be most useful 
in all Casimir problems involving nondisperssive media. In practical calculations the only formula needed is 



£y = (2--i)c(-«), 



1=0 



from which J2i=i v = — 1- Moreover, since the two last sums in (3.7) are known, 



" 1 1 o A 1 1 



E 

i=i 



72 = r 



E^I = 6 7r 



1=1 



we get, omitting the remainder in (3.7), 



E = 



64a 



1 



9 \ 



128 V2 



^ - 4 - 



423 (1 



16. 



16384 V 6 



^ - 16 



(3.8) 



(3.9) 



(3.10) 



The energy is positive, corresponding to a repulsive surface force. Remember, though, that we are working here with 
the nondispersive theory only. 



The structure of the three different terms in (3.10) is the following. The first term stems from the order 1/v in the 
Debye expansi on. T his term agrees with the results obtained by Milton |7| and Mil ton a nd Ng |sj in the case when 
the condition (2.13) holds. Numerically, the three terms between square brackets in ( ^To| ) are [1 + 0.06573 - 0.00610]. 
Thus the second term, stemming from the order 1/v 3 in the Debye expansion, describes a repulsive correction of 
about 6.6 per cent. Finally the third term, stemming from the order l/i/ 5 in the Debye expansion, describes a 0.6 
per cent attractive correction. We have thus improved the calculations in 0J^] by four orders in magnitude. The next 
correction, not included here, is of order 1/v 7 in the Debye expansion. 

Strictly speaking our consideration is applicable only when the condition eipi — 1 is satisfied. However keeping 
in mind a smooth dependence of the Casimir energy on the parameters specifying the field configuration, one can 
expect that the final formula can also be used in the case of a nonmagnetic dilute dielectric ball placed in vacuum (or 
a spherical cavity in an infinite dielectric surrounding) at least for estimation only. Therefore we present here some 
numerical estimations relating to sonolumincscence. With a good accuracy one can put in dimensional units 



E ~ 3^ 2 /(64a). 



(3.11) 



Take, as an example, |£| = 0.1, a = 4 • 10~ 4 cm. Then E ~ 2- 10~ 5 eV. This is immensely sm aller than the amount 
of energy (~ 10 MeV) emitted in a sonoluminescent flash. Furthermore, the Casimir energy (3.11), being positive, 
increases when the radius of the ball decreases. The latter eliminates completely the possibility of explaining, via 
the Casimir effect, the sonoluminescence for bubbles in a liquid. As known J26| emission of light takes place at the 
end of the bubble collapse. Recently an important experimental studies have been done to measure the duration of 
the sonoluminescence flash [^?J. In view of all this it is difficult to imagine that the Casimir effect, at least in its 
nondispersive version, should be important for the sonoluminescence phenomenon. 



IV. CONCLUSION 



Our method for calculating the Casimir energy E by means of the contour integral (2.7) proves to be very convenient 
and effective. As known, there are in principle at least two different methods for calculating E: one can follow a 



5 



local approach, implying use of the Green's function to find the energy density (or the surface force density). Or, one 



can sum the eigenfrequencies directly. Equation ( 2/? ) thus means that we have adopted the latter method here. The 
Cauchy integral formula turns out to be most useful in other context also, such as in the calculation of the Casimir 
energy for a piece wise uniform relativistic string A survey on this subject can be found in [E9]. The great 

advantage of the method is that the multiplicity of zeros in the dispersion function is automatically taken care of, i.e., 
one does not have to plug in the degeneracy in the formalism by hand. 

A remarkable feature of our approach is also the ultimate formula for the Casimir energy having the form of the 
spectral representation, i.e., of an integral with respect to frequency between the limits (0,oo) of a smooth function, 
spectral density. Evidently, for physical applications one needs to know the frequency range which gives the mai n 



contribution into the spectral density. An example of this representation for the partial energies Ei is Eq. ( 2.17 ), 
where the substitution y — ua should be done. As shown above, the partial energies Ei decrease rapidly as I 
increases. Therefore the most interesting is a few first values of I. In this case, as one could expect, the spectral 
density is different from zero when ua ~ 1. Keeping in mind the search for the origin of the sonolumincscence we 
put ||,|2(| a = 4 • 10~ 4 cm. Then the wave length of the photon in question turns out to be 25.0 • 10~ 4 cm, i.e., this 
radiation belongs to infrared region, while in experiments on sonoluminescence the blue light is observed fl26| . This 
fact also testifies against the possibility of explaining the sonoluminescence by the Casimir effect. 

It is worth noting that the spectral distribution of the Casimir energy is practically not discussed in literature while 



the space density of this energy has been investigated in detail (see, for example 18, (1983)]. From the physical point 
of view the space density and spectral density of energy in this problem should be treated on the same footing. One 
can remind here the treatment of the Casimir effect as a manifestation of the fluctuations of the vacuum fields |30| , 
these fluctuations being occurred in space and time simultaneously. 

It should be emphasized that in this paper we have neglected the dispersion effects when calculating the Casimir 
energy. Importance of this point has been demonstrated in JsT) . As for the elucidation of the sonoluminescence origin, 
we have to stress once more that in our consideration we have contented ourselves with the static Casimir effect only. 
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